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Abstract 

We present the threshold resummation for the cross section of heavy quark production 
in polarized deep-inelastic scattering to next-to-leading logarithmic accuracy and in single- 
particle inclusive kinematics. We expand our resummed result in Og to next-to-leading 
and next-to-next-to-leading order and study the impact of these higher order corrections 
on the charm structure function in the kinematical range accessible to the HERMES 
and COMPASS experiments. 
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1 Introduction 



Spin dependent deep-inelastic lepton-hadron scattering has received much interest during the 
past years, see Refs.|l|, ^ ^ Q for reviews. Presently, next-to- leading order analyses of the 
experimental data for the structure function gi allow for the extraction of polarized parton 
distributions. However, these are considerably less constrained than their unpolarized coun- 
terparts. This is particularly true for the gluon spin distribution A0g/p, which is currently 
determined from scaling violations. Numerous polarized hard scattering processes are sensitive 
to the gluon spin distribution, and the first moment of A(f)g/p contributes to the sum rule for 
the nucleon's spin content |[. 

Deep-inelastic leptoproduction of heavy quarks proceeds predominantly through the boson- 
gluon fusion process, and as such yields important information about the gluon distribution 
inside the proton. Polarized leptoproduction of charm quarks is therefore promising for an 
experimental determination of A(f)g/p. Presently, this measurement can be realized by the 
HERMES experiment [0. The upcoming COMPASS experiment is dedicated to measure 
A(f)g/p directly through the boson-gluon fusion process and subsequent polarized open charm 
production. COMPASS will probe both the region of photoproduction and deep-inelastic scat- 
tering and extend the kinematical range accessible to the HERMES experiment. 

In general, fixed target experiments, such as HERMES or COMPASS with centre-of-mass 
energies of ^/S = 7.4GeV and 14.1GeV respectively, operate rather close to the threshold 
for charm pair-production. In this region of phase space, the perturbatively calculable hard 
scattering cross section has potentially large higher order double-logarithmic threshold correc- 
tions. Thus, a resummation of these Sudakov logarithms to all orders in perturbation theory 
is required in order to regain control over the perturbative expansion. 

The purpose of this letter is to carry out this Sudakov resummation, which is done in sec- 
tion^. We work to next-to- leading logarithmic (NLL) accuracy (see Ref.0 for a recent review), 
and our final result provides an approximate sum of the complete perturbative expansion if, at 
each order, the Sudakov corrections dominate. Previously, studies of threshold resummation 
for polarized scattering processes have been performed for the polarized Drell-Yan cross section 
in Ref . [|To[] . By re-expanding our resummed cross section in ag to next-to- leading (NLO) and 
next-to-next-to- leading order (NNLO), we provide in section ^estimates of the size of the exact 
higher order corrections^ the calculation of which has not yet been completed fTl], |12|. Our 



estimates are valid in the region dominated by threshold logarithms and we study their impact 
on the charm structure function g^. In the kinematical region relevant to HERMES and COM- 
PASS, we investigate the sensitivity to the gluon spin distribution A(j)gjp and the dependence 
on the factorization scale. Section | contains our conclusions. 



2 Resummed differential cross section 

We begin with the definition of the exact single-particle inclusive (IPI) kinematics. We study 
the reaction of lepton (1) - proton (P) scattering, 

1(/) + P(j9) ^ 1(/ -q) + Q(pi) + XmV2) , (1) 

where Q is a heavy quark and X[Q] denotes any allowed hadronic final state containing at least 
the heavy antiquark. Neglecting electroweak radiative corrections]^, for longitudinally polarized 

^For the exact NLO corrections to polarized heavy quark photoproduction cf. Refs.[^ ^ |l5[ |l6|. 
^ QED radiative corrections have been studied in the leading-logarithmic approximation in Ref. |l7|. 
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protons, the difference of cross sections for reaction (|l]) may after summation over X be written 

TT|, ig, |T§ g 



as 



\2) 



dx dy d(j) dTi dUi dxdydcj) dTi dUi 



\2 ~ y — 



Q2 J dTidUi g2 dTidUl 



where a and M are the fine structure constant and the proton mass, and denotes the angle 
between the proton spin and the momentum of the scattered lepton, (Z — q) in Eq.(p. The 
lower arrow on d^a indicates the polarization of the incoming lepton in the direction of its 
momentum /, whereas the upper arrow on d^a represents the polarization of the proton, which 
is parallel or anti-parallel to the polarization of the incoming lepton. 

The functions d'^g^/ dTidUi, k = 1,2 are the double-differential deep-inelastic polarized 
heavy quark structure functions. The kinematic variables Q'^,x,y are defined by 

Q = ~q >0 , x = , y = . (3) 

2p ■ q p ■ I 

Additionally, we also have the overall invariants 

S = (p+qy = S'-Q\ T, = {p-p,f-m\ U, = {q-prY-m\ = S+T, + Ur+Q\ {A) 

Only the structure function gi is relevant asymptotically for M/Q ^ 1. It contains leading 
twist two operators and has a parton model interpretation. We therefore restrict ourselves to 
gi only. Furthermore, we neglect the contributions from initial state quarks and antiquarks. 
This is justified as they vanish at leading order. We consider only the gluon-initiated partonic 
subprocess 

l\q)+9{k) ^ Q(pi)+X'[Q](p'2), (5) 

where k = zp. The partonic invariants are 

s = {k + q^ = s' + q"^, U = {k - piY - , Ui = {q - piY - , Si = {p2f-m^. (6) 

The invariant S4 parametrizes the inelasticity of the partonic reaction in IPI kinematics. 
This distance from threshold is conveniently measured in terms of a dimensionless weight func- 



tion ws |21[ 

S4 2p2 - ks 2C ■ ks 

ws = — ^ — = • iV 

m 

The IPI kinematics vector ( is defined by splitting the recoil momentum = P2 + ks into 
the momentum p2 at threshold, and the momentum ks of any additional (soft) radiation above 
threshold. 

In double-differential form, d'^gi/ dTidU\ factorizes as 

2C?2^l(a;,Q^m2,Tl,^7l) } dz f 2^ a . ^2 2 2 / n\ /on 

S " y — A0g/p(z,/x ) Au;g(^-,S4,ti,Mi,(5 ,m ,/i . (8) 

The polarized gluon distribution in the proton is denoted by (^(^gjp and z is the momentum 
fraction with z~ = —Ui/{S'+Ti). The dimensionless function Aug describes the hard scattering 
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process and depends on the partonic invariants ti, ui and S4, defined in Eq.(H). The factorization 
scale is denoted by and throughout this paper taken equal to the renormalization scale. 

Sufficiently close to threshold, gi in Eq.(||) is dominated by those higher order contributions 
that contain plus-distributions of the type 



ln'(s4/m^ 
S4 




(9) 



They result from imperfect cancellations between soft and virtual contributions to the cross 



section. At order 0{a 



i+2\ 



0, 1, . . . we refer to the corrections with / = 2i + 1 as leading 



logarithmic (LL), to the ones with I = 2i as next-to-leading logarithmic (NLL), etc. 

The Sudakov resummation of these singular functions in Eq. @ uses the methods and results 
of Refs. pi], |2^ , p3| , |2^ . It rests upon the factorization of AcUg into separate functions Ai'g/g, 
S, and Hg for the jetlike-, soft, and off-shell quanta. This factorization, valid in the threshold 



region of phase space, implies a decomposition ||2^ of the total weight function w 



W = W^ + Ws 



Wi 



-Ml 



+ 



(10) 



which expresses the overall weight in terms of the individual weights and ws- Each of the 
functions Aipg/g{wi), S{ws), and Hg, which are conveniently computed in ( ■ A = gauge, 
organizes contributions corresponding to a particular set of quanta and thereby only depends 
on its own individual weight function. 

It is convenient to consider these functions in moment space, deffned by the Laplace trans- 
form with respect to the overall weight w, 



f{N)= ldwe-''^f{w) 



(11) 



Under the assumption of factorization the partonic cross section Acu^ in moment space can 
then be written in a factorized form, up to corrections of order 1/A^, as PTl ESI EH] 



/^Ug{N, ti, ui, fj,^, 



(12) 



i^'g(C,Q^m^a,(/i)) 



Ai,g/g{N^,k-C/fJ.) 

A^g/g{N^,fi) 



s 



m 

ml 



where = N{—ui/m^). The A^-dependence in each of the functions of Eq.(|T^) exponentiates. 
The jet and soft functions, Aipgjg and S in Eq. (p!^), can each be represented as operator 

. Dependence on the spin degrees of freedom can be kept 



matrix elements p2|, ^ 
explicit in these elements, such that the methods of Refs.[^, 22, 23, 2|] for the resummation 
of the singular functions in Eq. via appropriate evolution equations can be generalized in a 
straightforward way to account for initial state polarizations. The jet-function Aipg/g obeys two 
evolution equations, namely the renormalization group equation and an equation describing the 
energy dependence via gauge- dependence |^ |2^, Solving both equations sums Sudakov 
double logarithms. Mass factorization introduces the density Acpg/g as a counterterm and 
requires the choice of a scheme. In the MS-scheme A(j)g/g has no double logarithms. The 



soft function S obeys a renormalization group equation, which resums all NLL logarithms 
originating from soft gluons not already accounted for in Aipg/g. In general, the soft and 
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the hard function S and Hg depend also on the colour structure of the underlying scattering 
reaction, but for deep- inelastic production of heavy quarks this dependence is trivial [^. 

The final result for the hard scattering function Aug in the MS-scheme in moment space 
resums all In in single-particle inclusive kinematics. We obtain 



AcD^*^^(iV,^l,Ul,Q^m^/i^a,(/i)) = 

i^,(C,Q^m^a,(m))5(l,C)exp|2 



(13) 



X exp 



dw 



w 



(1- 



-Nuw 



^A<^iasiVX2k ■ 0) + ]-i^giasiw2k ■ ()) 
A 2 



X 



m/N 2k-(^ 

expj J ^2Rer5(«.(/i'))|exp|-2 | ^ (7^ (a,(^')) - 7./.(iV., «,(/i'))) 



where is defined below Eq.(|^). To NLL accuracy, as defined below Eq.(^, the product 
HgS on the second line of Eq.(^) is determined from matching to the Born result at the scale 
H = m. To this accuracy the product HgS is also insensitive to the choice of treatment of 75. 

The second exponent in Eq.(|l^) gives the leading iV-dependence of the ratio Aifjg/g/ Acpg/g 
with 



Ag{as) 



a 



Ca- + -CaK[^\ +... 



a. 



TT 



a. 



i/S(a,) = 2Ca— + ..., 



IT 



where K = (67/18 — n^/G) — 5/9nf can be inferred from Refs. P7|, using Ref. 
Ref. [0). The scale evolution of the ratio Aipg/g/ A(j)g/g is governed by 

a. 



79 ("s) 

-fg/giN, as) 



62— + 

vr 
a. 



TT 



a 



-{CAlnN -b2)-[ — ] [-CAK\nN] + ... , 



TT 



with "fg/g calculated in Refs. ||27|, |28| , and the soft anomalous dimension to order as is 



as 



TT 



Ca 



Ca 



-{[^-CF]{Lp + l)-^\\n 



+ ln- 



m 



+ . . . 



(14) 
(cf. also 

(15) 
(16) 

, (17) 



with/?= ^1 -4mVsandL/3 = (l-2mVs){ln(l-/5)/(l+/9)+i7r}//?. We note that the various 
anomalous dimensions in Eqs. (|I^)-(|T^) are identical to NLL accuracy to the corresponding 



quantities of soft gluon resummation for unpolarized scattering |25 



Eq. ([13D represents the central result of this paper. It provides the sum of Sudakov logarithms 
due to soft gluon emission to all orders in the perturbative expansion and accurate to the next- 
to-leading logarithm. 



3 Finite order results 

We expand our resummed result in as up to second order so as to provide NLO and NNLO 
approximations to NLL accurracy for the partonic single-particle inclusive double-differential 
cross section difference Aa^*g. 
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For the process 7* + — > Q + Q, the single-particle inclusive double-differential cross 
section difference is related to the hard function Aug of Eq.(^). To NLL accuracy it can be 
written in a factorized form as 



,2 C?^Acr^*g(s , ti, Ml) o 2 \ f . ,n2 2 2 /■ \ 



dti dui 

jBorn / 



(18) 



with the Born level hard part Ai?B.°™ given by [^O], ^ 



A q 



- + — + — + — +2gM- + _ + 

Ml ti \ti Ui Uf tf J \ti Ui S' 



The NLO and NNLO soft gluons corrections to NLL accuracy are 



ln(s4/m^) 



2C, 



S4 



+ 5(s4)CAln 



-Ml 





■ 1 ■ 




+ 










+ 








) In 









C^(ln(^)+ReL^-ln(i^ 

.Ml/ 



(19) 



(20) 



2C^(ReL^ + 1) 



"In^(s4/m2) 



S4 



+ 



In^(s4/m2) 



Si 



3Cl In 



ln(s4/m 
S4 



(^) 



ReL/3 — In | — - 



In ( ^ M Cl ( -2 In 



-2C^(62 + 3CF(ReL^ + l))| + 
-2ReL^ + 2 In (^^) + In f^"! + 2 (62 + 2 (ReL^ + 1)) 



+ 



1 

S4 



,2 / 



lnM^H-Clln(-;;i) -=:CAh 



-Ml 



(21) 



with fi the MS-mass factorization scale, 62 = (HCa — 2?2j)/12 and given below Eg. (P!7D . 
We have checked that to NLL accuracy, the result in Eg. pO|) agrees with the exact 0{as) 
corrections calculated in Refs. ITII, 0^ . 



Next, we perform a numerical investigation of the results (|20| ) and (pi]). To that end, we 
relate the total partonic cross section difference to dimensionless coefficient functions Ac^^''-*, 



a tts e, 



2 00 



k=0 



k 

— "V 

1=0 



^E(4vr«^(/^))'EAf'')(r/,01n' 



I 



m 



2 ■ 



(22) 



As completely inclusive guantities the functions Ac^l^''-* depend only on the scaling variables 




10"^ 10"^ 10"' 1 10 10^ 10"^ 10"^ 10"' 1 10 10^ 



ri = s/(4m^) - 1 ri = s/(4m^) - 1 

Figure 1: (a) The ?7-dependence of the coefficient functions Acg^'°''(r/, A; = 0, 1 for = lOGeV^ 
with m = 1.5 GeV. Plotted are the exact result for Acg'^'^^ (solid line) and the LL approximation 
(dotted line) and NLL approximation to Acg''^^ (dashed line), (b) The ry-dependence of the coefficient 
function Ac^'^\r],^) for = lOGeV^ with m = 1.5 GeV. Plotted are the LL approximation (dotted 
line) and NLL approximation. 



where rj is a measure of the distance to the partonic threshold. 

In Fig.[^ we plot all coefficient functions Ac^'^'°);fc = 0,1,2, i.e. those not accompanied 
by scale-dependent logarithms for = lOGeV^ and m = 1.5 GeV. Only the Born function 
/\c^^fi) in Fig.|I]a is known exactly For Ac^^'^^ in Fig.^a and Ac^^'°) in Fig.[^ we can 

give estimates to LL and NLL accuracy. In the LL case, we keep only the [In(s4/m^)/s4]+ and 
[In^(s4/m^)/s4]+ terms in Eqs. (pO]) and (|2T| ) respectively. We find both Ac^^'°) and Acf'^^ to be 
sizable near threshold where the large logarithms dominate, while they tend to be numerically 
rather small at larger values of rj. Moreover, we also investigated some next-to-next-to leading 
(NNLL) logarithmic contributions^, such as the Coulomb corrections |Tl|, ^ for Ac^^'°) 
and for Ac^^'°) those NNLL terms, which we obtain from the expansion of the resummed result 
Eq.(13). Numerically, we found these NNLL terms to have an effect on Ac^^'"^ or Ac^^'"^ of the 
order of 5% as compared to our NLL corrections. 

Since no exact results for Ac^^'°) and Ac^^'°^ are available we are unable to judge the quality 
of our approximation. However, similar investigations in the case of unpolarized heavy quark 



production p5[, where exact one-loop results are known suggest that generally NLL accu- 



racy provides a very good approximation of the true result for not too large^, because K^^\ 
are the same as for the unpolarized case, only differs from AB^?™. 

Let us now turn to those coefficient functions multiplying scale-dependent logarithms. These 
are Ac^g'^^ in Fig.|]a and Ac^g'^\ Ac^^'^-' in Fig.0b, which we plot for the same values of para- 
meters as in Fig.|l|. We find in particular the estimates for Ac^^'^^ and Ac^^'^^ to NLL accuracy 
to be large near threshold. Additionally, we are able to obtain the exact result for the function 

^ For complete resummation to NNLL accuracy one would need to match the resummed result Eq .(|l^) a t 
NLO, which requires knowledge of all one-loop soft and virtual corrections. These are not yet available |ll|, |l2|. 

^In the regime of large Q^/m^ on the other hand, the coefficient functions Acg'^'"-' can be approximated with 
different methods based on the operator product expansion |20|. 
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Figure 2: (a) The r/-dependence of the coefficient function Acg^'^\r],^) for = lOGeV^ with 
m = 1.5 GeV. Plotted are the exact result (solid line), the LL approximation (dotted line) and the NLL 

(2 I) 

approximation (dashed line), (b) The ry-dependence of the coefficient functions Acg (??,^), / = 1,2 
for = lOGeV^ with m = 1.5 GeV. Plotted are the LL approximation (dotted line) and the NLL 
approximation (dashed line). 




Ac^^'^^ by means of renormalization group methods. 



1 

Ac(^'i) {v{x),0 = dz {h,5{l-z) - ^APW(^)) Acf°) [nQ.dj , (24) 



with a = 1 + Am^ IQ"^ and the one-loop splitting function ^Pg^g given in Refs. 3^, 
Fig.|]a clearly shows that the approximation based on NLL accuracy traces the exact result 
(|2^) extremely well even at larger r]. For Ac^^'^^ we have checked that NNLL terms obtained 
from the expansion of the resummed result Eq.(|13|) have only a small effect. 
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Figure 3: (a) The polarized gluon distribution rjA(j)g^p{7],x, fj.'^) for x = 0.1 and /j = 1.5 GeV. 
Plotted are the parametrizations GS A of Ref.|§ (sohd), GS C of Ref.|3| (dashed) and GRSV 
valence of Ref.p^ (dotted), (b) Same as (a) for x = 0.01. 



Let us now investigate the impact of the approximate higher order perturbative corrections 
on the inclusive hadronic structure function gi. In terms of coefficient functions, the charm 
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Figure 4: (a) The x-dependence of the charm structure function xg\{x,Q'^ ,m?) with the gluon 
distribution GS A of Ref.@ for /i = m = 1.5 GeV and = lOGeV^. Plotted are the results at 
leading order (solid), at NLO to LL accuracy (dotted), at NLO to NLL accuracy (dashed), at NNLO 
to LL accuracy (spaced dotted) and at NNLO to NLL accuracy (spaced dashed), (b) Same as (a) 
with the gluon distribution GS C of Ref.|36|. 



structure function gl can be expanded as follows: 

gl{x,Q\m^) (25) 

on III X J in 



ax 
A 



2n 



2 



f d{\og^ T,) In 10 r/A0,/p(r/, x, /i^) E(4vr«. (/i))' E ^4'''^ In' ^ 



-oo 



A;=0 Z=0 



where a = l + 4m^/Q^, A = log^o(^{V^~ — 1), and ^(pg/p represents the polarized gluon 
distribution. We recall that the MS-scheme has been chosen and that contributions from light 
initial state quarks are neglected. For A0g/p we compare the parametrizations of Refs. [p^, |37 



For our analysis at NLO (LO) we use a 2-loop (1-loop) running coupling with rij = 3 light 
flavours, a charm (pole) mass of m = 1.5 GeV ^2[, and Aqcd = 0.232 GeV. 

In Eq. (^5|) , we have chosen to rewrite the standard expression for the convolution of g\ in 
a form, that facilliates the investigation of partonic threshold effects on g\. To do so, we plot 
r]A(l)g/p as a function of t] in Fig.§, and examine the support it gives to the coefficient functions 
in the second expression of Eq.(p5D. By comparing Fig.|^ with Figs.|l] and ^ one can directly 
judge over which ranges of t] the function 't]A(j)g/p becomes large and hence samples the partonic 
coefficient functions. 

Fig.|^ reveals that for all parametrizations this happens indeed in the threshold region for 
x>0.01 and at scales around ~ lOGeV^. We therefore expect our estimates for Ac^^'^^ 
to provide a good description of the true higher order corrections for gl for x>0.01. The 



gluon densities from set A of Ref . and the valence scenario of Ref . are both positive and 



similar in shape, while the gluon density from set C of Ref. |^ relaxes the positivity constraint 
on A(f)g/p and oscillates.0 



^Recent results from the HERMES collaboration |38 indicate a positive ratio of polarized over unpolarized 
gluon distribution A(j)g/p/(j)g/p at intermediate x. 
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)x[GeV] |x[GeV] 

Figure 5: (a) The factorization scale dependence of the charm structure function 

gi{x,Q'^ .m^ , li'^) / gi{x,Q'^ , ji^ = vn?) with the gluon distribution GS A of Ref . |3|] for m = 1 . 5 GeV , 
= lOGeV^ and x = 0.05. Plotted are the results at leading order (solid), at NLO to LL accuracy 
(dotted), at NLO to NLL accuracy (dashed), at NNLO to LL accuracy (spaced dotted) and at NNLO 
to NLL accuracy (spaced dashed), (b) Same as (a) with the gluon distribution GS C of Ref. p^] . 



Next, we investigate the x-dependence of g'l. In particular we are interested in a comparision 
of our NLO and NNLO estimates (|20| ) and ( pT] ) with the leading order result |3y, In 



Fig.^a we plot the leading order result for xgl as well as the NLO and the NNLO approxima- 
tions over a range 0.007 < x < 1. We choose = 10 GeV^, a fixed value of the factorization 



scale fi = m = 1.5 GeV and the gluon parametrization GS A of Ref.[36|. We obtain similar 



results with the gluon distribution in the valence scenario of Ref. . We find that the pertur- 
bative corrections are sizable, both in the region of intermediate x, around 0.05 and at smaller 
X, where however, the approximation is less certain to work well. To assess the quality of our 
approximation, we compare at each order our results to LL and to NLL accuracy. At interme- 
diate X, the small differences between LL and NLL accuracy show a very good stability of the 
threshold approximation for the description of g^. Towards smaller x however, these deviations 
increase and indicate that our approximate methods fail for x < 0.01. 

In Fig.^b we repeat the plot of xgl for the same parameters as in Fig.^ but with the gluon 
density GS C of Ref. . Again, the perturbative corrections are sizable over the whole range 
in X, but the shape of the curves is completely different. Both the coefficient functions and the 
gluon density oscillate, which leads in particular at intermediate x around 0.1 to destructive 
interference, with gl being only marginally different from zero. 

The final issue we study is the dependence of gl on the mass factorization scale. In general, 
leading order calculations exhibit a strong sensitivity on the factorization scale, which is usually 
reduced once higher order corrections are taken into account. In addition, there are general 
arguments supporting a reduction in scale dependence from including soft gluon effects |]40[1 . 
Therefore, we are interested in the effect of our NLL-approximate NLO and NNLO results (|20| ) 
and (|2T|) on the scale dependence of g^ in comparision to the leading order calculation [^, 
Note however, that the arguments of Ref.[^] leading to genuine NLO soft gluon 
resummations rely on NNLL accuracy, which is not yet available for g^, see footnote ^. 

We plot the ratio gl (fi) /gUfi = m) over a range m < fi < lOGeV and fix m = 1.5 GeV, 

= 10 GeV^ and x = 0.05. In Fig.^a we use the parametrization GS A of Ref. [p6l , but the 
following conclusions hold also for the gluon density in the valence scenario of Ref.[^. We 
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find at NLO that the corrections based on LL accuracy only are not sufficient to reduce the 
scale dependence. Clearly, soft gluon effects need to be approximated at least to NLL accuracy 
to achieve the desired result, a feature that has also been noticed in studies of unpolarized 
heavy quark production Our best estimate for makes use of the exact expression for 
the coefficient function Ac(^'^) of Eq.(|4D. However, since it differs only slightly from the NLO 
result to NLL accuracy for the chosen parameters in Fig.^, we do not display it here. 

In Fig.^ we repeat the analysis for the parametrization GS C of Ref.[36|. In this case. 



the NLO and NNLO approximations do not reduce the scale dependence. This is due to the 
oscillating gluon density, which leads to gl being close to zero at the x-value chosen, and even 
causes g^ to change sign, depending on the scale. 

To summarize, the study of the charm structure function as in Figs.^ and |^ shows that 
the soft gluon estimates of higher order corrections to the coefficient functions are well under 
control and give stable predictions for gl at scales not too large and x>0.01. On the 
other hand, in the chosen kinematical range the effects of higher orders do not upset the 
sensitivity to the polarized gluon distribution function, with different gluon parametrizations 
leading to qualitatively different behaviour for gl. Therefore, measurements of the charm 
structure function allow to further constrain the polarized gluon density. Finally, we found 
that the sizable variation of gl at leading order due to different values of the charm mass are 
not reduced by our approximate higher order corrections. 



4 Conclusions 

In this letter we have investigated the effects of soft gluons in the case of polarized lepto- 
production of heavy quarks. Our study helps to keep theoretical control to next-to-leading 
logarithmic accuracy over large higher order corrections due to threshold logarithms and we 
have presented analytical results for the all-order resummed cross section in single-particle 
inclusive kinematics. 

Moreover, we have provided both analytical and numerical results for the approximate NLO 
and NNLO perturbative corrections to the deep-inelastic charm structure function gl, to NLL 
accuracy. As a consequence of the scattering reaction being largely driven by initial state gluons 
with not much more energy than required to produce the final state, we have found our analysis 
of gl to be well applicable in the kinematical range accessible to the HERMES and COMPASS 
experiments. In this region, our NLO and NNLO estimates can help to reduce theoretical 
uncertainties and may assist in the theoretical interpretation of future gl measurements. 
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